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Abstract 

From the minimal action principle follow the Hamilton equations of evolution for geomet- 
ric optical rays in anisotropic media. As in classical mechanics of velocity-dependent poten- 
tials, the velocity and the canonical momentum are not parallel, but differ by an anisotropy 
vector potential, similar to that of linear electromagnetism. Descartes’ well known diagram 
for refraction is generalized and a factorization theorem holds for interfaces between two 
anisotropic media. 


1 Fermat’s principle 


Fermat’s principle states that the light ray joining two points in an optical medium takes the path 
where it employs an extremal time [1]: 

6 f dt = 6 f dsn(q(s),q(s)) = 0. 

Ja Ja 

Here we denote by ds the length element along the ray q, the ray direction by q = — - and by n the 

ds 

refractive index of the medium. The refractive index characterizes the optical medium. Constant 
n indicates that the medium is homogeneous (invariant under translations) and isotropic (invariant 
under rotations). In anisotropic media, the refractive index depends also on the direction of the 
ray [2]. 

We use one of the Cartesian coordinates of 7 Z 3 as the evolution parameter to describe the 


evolution of the ray q = 
Fermat’s principle as [3] 


q 

z 


_ . da dz 

Defining v = — with ds — . . 

dz v 1 — q 2 


dz \/l -f v 2 we write 


rzB 

8 dzL(q(z),z]v(z)) = 0, 

Jz A 


with the Lagrangian function L(q, z;v) = \/l + v 2 n(q, z\ v). 
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2 Evolution equations 

The Euler-Lagrange equations that follow from the Fermat principle are [4] 

d_ _dL 
dz dq ’ 

where the canonical momentuffi is 

dL nv , /n — ■ , m ■ \ 
P = , f + VI + v 2 — = nq + A(q,z,q), 

d\ V 1 + V 2 crv 


and we define the anisotropy vector 


A = \/l + v 2 


dn 

fa 




We obtain the Hamilton evolution equations through the Legendre transformation 


dq _ dH dp_ dH_ 
dz dp ’ dz dq ’ 


with the Hamiltonian function 

H(q,z] p) = p • v - L(q,z\ v) 


- yjn 2 - |p - A| 2 + 


(p - A) • A 

- |p - aF ' 


In anisotropic media, the three-vectors of ray direction q, momentum p, and anisotropy A , are 
thus characterized by: 

p= nq + A(q,q), p z = -H, \q\ = 1, \p - A\ = n{q,q), 


i.e., we have the orthogonal decomposition of momentum p into ray direction q and the anisotropy 
three- vector 


-+ 

A = 


V’-n 


bl=i 




(i -?r) 


dq 


141=1 


A dn 

e, w + 


dn 

sin 6 d<f) 


The anisotropy vector is orthogonal to the direction of ray propagation q. 

While |(f | sweeps over the ray direction sphere S 2 , the vector p — A draws out a closed surface 
n(q, q) — the ray surface at the space point q, and the three- vector p ranges correspondingly over 
another closed surface that we call the Descartes ovoid of the anisotropic medium at q. 

The Hamilton equations are thus written in manifestly euclidean-covariant form as 


dq dH p — A dp dH _ n dn 

dz dp p z — A z ' dz dq p z — A z dq 

dn -* 

From the second equation it follows that dp x — = 0. As in the isotropic case, we get the Ibn 

dq 

Sahl [5] Snell law of refraction between two anisotropic media [6]. 
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FIGURE 1. Dipole medium: the momen- 
tum p of a ray is obtained by adding the 
direction vector q times n° to the dipole 
vector of the medium D. The ^nisotropy 
vector A ranges over a cardioid-type sur- 
face. 



3 Dipole anisotropic media 

Consider the refractive index with linear dependence on ray direction 

"(9,9) = »°(9) + D(q,q), 

»(?.?) = £ D,(q)'q, = D(q) T 'q. 

j=x,y,z 

We call n° the monopole part of the medium and D its dipole vector. The anisotropy vector is 

A (l) = (1 - qq T )D = D - D T (q)q = (q X D) X <f. 

This vector lies in the plane of q and D, and is orthogonal to the ray direction q. The relation 
between ray direction and optical momentum is p — n°q + D . While q £ the Descartes ovoid 
is a sphere of radius n°(q) and center at D . (See Figure 1). 


4 Quadrupole media 

Consider now a refractive index with quadratic dependence on ray direction q 

n(q,q) = n°(q) + Q(q,q), 

W, 9 ) = J2 Qi,k = q T Qq- 

j,k= X} y,z 

We have a ^-quadratic summand with coefficients Qj^ in a 3 x 3 optical quadrupole matrix Q 
that must be symmetric and traceless. It is common to restrict consideration to principal axes ; in 
that frame of reference, Q = diag (Q xy Q y , Q z ) and Q x + Q y + Q z = 0. Then, the anisotropy and 
momentum vectors are 

A^^2(\-'q'f)Qq = 2[Qq-Q{l'q)'ql 

P= [n(q,q) + 2(1 - qq T )Q]q = (n° + 2 Q- q T Qq)q. 

In two-dimensional optics, 2x2 symmetric traceless matrices have two independent coefficients 
that describe the ellipticity and orientation of the figure. When q ranges over the sphere of 
directions «S 2 , p will range over the Descartes ovoid of the quadrupole medium. (See Figure 2). 
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FIGURE 2. n°q ranges over the sphere 
(circle in the two dimensions of the figure) 
of radius n°, nq over the peanut-shaped 
surface, and the momentum p draws a 
Descartes oval. The thin lines joining 
points on the circle and on the oval re- 
late the direction of the ray with the corre- 
sponding direction of the momentum vec- 
tor. 



5 Free propagation in homogeneous uniaxial media 

For the uniaxial quadrupole media we can write the refractive index as 

v 0 

n(q) = n° + (x, y, z) | 0 v 0 

0 0 -2v J 



where n° is the monopole part and v is a quadrupole anisotropy coefficient. 
Putting q in terms of v, we can write out the components of momentum as 


p = (n° + 4 1 /) q - 3i/q 2 q = (1 + V 2 ) 3/2 [(n° + Au) + (n° + ^)v 2 ]v 


For free propagation, the Hamilton equations and their solutions are: 

J = V => q {z) = q(0) + 2 V, ^ = 0 =* p (z) = p(0). 

Although the solutions are apparently independent of the anisotropy of the medium (they are 
straight lines in space), the anisotropy is expressed through the relation between the ray momen- 
tum p and the ray direction v. In isotropic media, the momentum vector is n times the direction 
vector and we can easily invert this particular case to [7] 


v = = = — , | v | = n° tan 0 (v = 0). 

^/(n 0 ) 2 - p 2 Vz 

In the more general uniaxial anisotropic case, to find a simple closed inversion, we expand this 
equation with a Taylor series in (v 2 ) fc v for k = 0, 1,2, . . . and propose a similar expansion of v in 
powers of (p 2 ) fc p. Equating the series we find the expansion coefficients 


v( P ) = 


P + 


-n 


° + 5i/ 


7 P P + 


|(n 0 ) 2 + yn°i/ + 51i/ 2 


n° + Au (n° + 4 u) 

^(n°) 3 + j-(n°) 2 t' + 114n°i' 2 + 650i/ 3 
+ ( n ° + 4 1/) 10 


(n° + Ai >Y 
(P 2 ) 3 P + 


(P 2 ) 2 P 
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FIGURE 3. Descartes diagram to construct 
refraction angles between two anisotropic 
media. 



The evolution Hamiltonian is then 

H( p) = p v — (n° — 2v)\/\ + v 2 — 3^ 


+ v 2 


1 2 n° + 101/ 2 2 3(n 0 ) 2 + 60nV + 408!/ 2 2 3 

= - 2 " ) + s? p + w (p) — — (p) 

, 5(n 0 ) 3 + 150(n°) 2 i/ + 1824nV + 10 400 v 3 ,_ 2yl , 

T ^ an / •MO IP 1 ' ' ‘ ■> 


128(n e ) 


where n e = n° + \v plays the role of an effective paraxial refractive index. In the isotropic case 
(when v = 0), this is the expansion of --\/n 2 — p 2 , the well-known optical Hamiltonian for such a 
media [8]. 


6 Finite refraction 


Let us consider the case when rays cross the flat interface z — 0 between two different aligned 
uniaxial quadrupole anisotropic media. Let their refractive indices in two half-spaces be n(q) 
and n'(q) with monopole parameters n° and n°\ and quadrupole anisotropy coefficients v and 
j/, respectively. The refraction law claims that the projection of the momentum vector on the 
refracting surface is conserved , which for our refracting surface gives p = p / . Generally, the 
incident and refracted rays are not coplanar with the surface normal. However, in the aligned 
uniaxial case both refractive indexes are axially symmetric (under rotations around z-axis) and 
the two anisotropic vectors are coplanar with the suface normal (z-axis). Refraction in our case 
is thus coplanar. 

Using ‘ruler, compass and plotter’ on the plane Figure 3, we construct the Descartes diagram 
for the point at the interface joining two ‘half Descartes diagrams’ and matching the length of the 
momentum vectors p and p' on the interface. To find the angle of refraction 9' in terms of the 
angle of incidence 9, we construct q'(q; n°, v\ n 0/ , v') expanded in series of sines, and find 


• n . n 3 

sin 9 — — 7 sin 9 4 ; 

n e n e 


n 


n c 


sin" 


_ 27 f n e \ 2 , fn e \ 3 , 


sin 


9 
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The first summarid is the very well know law of sines (Ibn Sahl-Snell law [5]); it is here also the 
paraxial approximation with the ratio of effective refractive coefficients. The succeding terms are 
corrections of orders u k and sin 2 * +1 0 due to anisotropy. 

7 The roots of refraction 

We consider now the ray transformation due to refraction at a smooth surface S(q) = £(q) — z = 0 
between two general anisotropic, homogeneous media n(q ) and n'(q). The rays in the first and 
in the second media are given correspondingly by the equations 

q(z) = q + zv, p(z) = p, z<z, 
q[(z) = q' + W, p \z) = p', z > z, 

where we have indicated the point of impact at the refracting surface by bars q — (q, z = ^(q)). 
We can formally consider the second pair of equations also on the left of the refracting surface, 
z < z. It allows to parametrize the rays behind the surface by the coordinate q' and momentum 
p' on the same screen z = 0; v and v' are the two ray directions on the screen. Thus, the point 
of impact coordinates can be written in two ways: 

q(z) = q + C(q)v = q = q' + C(q)v' = q!{z). 

This is the first root equation of refraction [9]; it is an implicit equation for q. 

The second root equation follows from the conservation of the tangential component of mo- 
mentum and implies the refraction law. If the normal to the surface S is denoted by VS(q) = 
((*, C».-l) = (S(q).-l) then we have (p — p f ) x V5 r (g) = 0. As we know, the momentum 
vector has components p = (px.Py^Pz) = (p, — H). Denoting the Hamiltonians before and after 
the refracting suface as H and H f we can rewrite the last equation containing the vector product 
as 

P - # (p)E(g) = P = P' - #'(P ')£($). 

This is the second root equation determaning explicitly p once q has been found. 

We have thus determined the root transformation for generic surfaces S — {(q) — z = 0 between 
homogeneous, anisotropic media. On optical phase space the root transformation is 

: q *-*■ q = q + v(p)C(q), 

: P ^ P = P-^(p)S(q), 

where v(p) and H( p) contain the refractive index function n(q). From our construction follows 
that the refracting surface transformation 

Sn,n';< = (q, P) («l',p / ) 

thus factorizes into the product of the root transformation in the first medium and the inverse 
root transformation in the second medium, <S n ,n';C = (^n' ; c) _1 - When the surface S' is a 
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FIGURE 4. Refraction at a surface is 
a map between phase space points (q,p) 
and (q , ,p / )- This transformation visi- 
bly factors into transformations back and 
forth from the point of impact q on the 
surface z = £(q). 



z = constant plane, the second root transformation is simple free flight by generic 2 . The root 
transformation is illustrated in figure 4. 

Let us consider explicitly the example of the symmetrical surface under rotations around z-axis 

C(q) = C 2 q 2 + C4q 4 + • ■ • • 

The refraction by such a surface is determined to third aberration order as [9] 

q' = 

P' = 


The paraxial part of the transformation is recognizably that of a quadratic surface. 


q - C 2 


1 


1 


n°> -f 4 u' n° + 4v 


2 , . „,2 1 "°' - 2,/ ! - K - 2 "] 2 
qp + 2(2 7FTT* qq ’ 


p + 2C2([n°' - 2 *'] - (n° - 2 ,])q + (2 (^57^; “ P 2q 


, 0 / 


-4C: 


2 L rr' — 2v f \ — [n° — 2v\ 

2 


n 0/ + \v 


p qq - 2C 2 


n 0/ — 2v'] — [n° — 2v] 
n 0> -f 4 v 


q p 


+4 ^3 ([„»' _ 2 v'\ - K - 2p])j q 2 q* 
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